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ABSTRACT 


In this thesis we consider the lot sizing of components 
in a Group Technologic (GT) production system* The manufac- 
ture of components in the GT production environment involves 
the establishment of cells comprising of several machines. 

Each cell is capable of handling a part family and the 
components move through the various machines in the call as 
far as possible sequentially. In the present work we assume 
that it is possible to obtain a part family and a production 
cell so that no back tracking of components is required. 

Further, it is assumed that each stage comprises of one machine 
A review of the literature suggested that in such a production 
environment, work-in-process (iflP) inventory contributes sig- 
nificantly to the total cost of production/inventory system. 

In the present work we have developed four mathematical models,- 
for lot sizing of components belonging to a part family which 
is to be manufactured in a GT cell. 

Kie mathematical models developed assume a constant 
demand rate for the components, and no inter-cell movements 
for all components. 

The first model basicall3)" is an extension of the Wilson's 
model for determining the economic production quantity. The 
various machines are identified as stages and a component 
belonging to the part family is manufactured sequentially on 
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the various stages. The model assumes that the lot (economic 
production quantity/-) will move from one stage to the next 
stage only after the entire lot has been maniifactured on the 
present stage. It is shown that ¥IP inventory influences the 
economic lot size as well as the total cost of the productiojq/ 
inventory system. 

The second model considers the splitting of the lot at a 
stage into batches to reduce the contribution of the ¥IP inven- 
tory. Each batch after it is manufactured at a given stage is 
transported to the next stage for further processing. The 
model is named as constant lot size model with lot splitting. 

A two stage algorithm involving a heuristic procedure and an 
optimisation procedure has been suggested for this model. 

The third model called the variable lot size model 
assumes the Crowston integrality Theorem, on the lot sizes at 
various stages, Crowston’ s integrality theorem states that 
if i denotes any stage, a(i) denotes its successor stage and 
N denotes the final stage then there exists a set of optimal 
lot sizes such that for all i < N the ratios 

= ^i/^a(i'' positive integers, A dynamic programming 

algorithm has been suggested for solving the variable lot size 
model. 

In the last model, we combine the importtmt features of 
second and third models, viz,, the splitting of the lot at a 



stage and -assigning variable lot sizes for the lots to be 
produced at various stages. This model is referred as the 
variable lot size model with lot splitting, A heuristic 
solution methodology has been suggested for this model. 

The solution methodology of the second model is ex- 
plained through an illustrative example. The proposed algo- 
rithms are coded in Fortran - 10 for implementation on 
DEC 1090 system. The coii^^utational performances of the 
various algorithms have been investigated for problems of 
varied sizes. 

Problem size varied from single stage to 10 stages. 
The input parameters for various problems were generated ran- 
domly in the specified ranges. The constant lot size model 
with lot splitting is compared with variable lot size model 
with lot splitting. 



CHAPTER I 


INTRODUCTION 

1.1 GROUP TECHNOLOGY CONCEPT; 

Group Technology (GT) is a very progressive method of 
organising production and especially it is becoming popular in 
those industries viiich are engaged in medium and small hatch 
production. The principles of mass production applied in 
hatch production would result in lot of duplication of paper 
work and loss of machining time. One fundamental step capa~ 
hie of bringing the mass production principles within the 
reach of mediimi hatch production is the adaptation of group 
machining approach. Professor Mitrafanov ( 1 ) is the pioneer 
in the field of group production. 

Group Technology is a technique of identifying and 
bringing together related or similar components in a production 
process in order to take advantage of their similarities in 
design, dimensions, geometrical shapes, raw material and 
tooling, by making use of the inherent economies of flow 
production methods. The aim is to substantially reduce the 
set'^up times and to improve the delivery performance by 
reducing the throughput times. This is achieved by organis- 
ing a large number of diverse components into families which 
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require similar manufacturing processes and providing the 
most suitable manufacturing facilities for groups of families 
by designing cells having machine centres exclusively for 
processing the group of components. The physical nature of 
difference betv/een traditional batch production and group 
production is shown in Fig, 1,1. 

The introduction of group concepts in production 
organisation has two main origins. Firstly methods were 
developed by engineers who were mainly interested in finding 
the techniques which would reduce the set-up time resulting 
in increased capacity and better utilisation of equipment. 
Secondly, behavioural scientists advocated the use of grovip 
concepts in organisational design for increasing the workers 
motivation and ;3ob satisfaction. The models based on group 
concepts developed by the engineers and behavioural scientists 
have number of similarities in terms of structure and solution 
methodologies and have been implemented successfully in varied 
situations. 

Burbidge ( 2 ) has listed number of successful applica- 
tions of group technology'’ in industry. I-iaximum benefits of 
GT applications result in a production environment involving 
a wide variety, medium volume of production of components, 
req’uiring different types of production facilities. 

The following benefits result from the introduction 


of Group Technology, 
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1* For each component all the required operations are 
done inside one cell. This results in the reduced throiaghput 
time, 

2, It is possible to sequence the loading operations well 
in advance since the components (parts) included in the 
family are prespecified. Thus it becomes possible to obtain 
a better load balance for the production cells and schediiling 
of components on machines within the cell, 

3 * By the application of composite component principle 
tooling, jigs and fixtxires can be standardised, there by 
reducing the set-up times considerably resulting in increased 
available capacity. 

4. Duplication in design and process planning effort can 
be reduced due to simplification in variety of design and use 
of schemes of component classification and coding. 

5. The responsibility for quality can be assigned as all 
the operations are performed in the cell. Hence quality levels 
improve for all components. 

6. GT simplifies the production control, material flow 
system and material handling. 

7. %)art from operational benefits, the social and psycho- 
logical benefits include better motivation, higher job satis- 
faction, improved skills and better quality of life. 
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1,2 DESIGN OF GT PRODUCTION SYSTH^: 

The design of GT production system embraces several 
functions such as component classification and coding, estab- 
lishment of groups and cells, development of a group layout 
and determination of lot sizes, sequences, schedules for com- 
ponents and load on cells. 

The following are the main characteristics of a GT 
production system. 

1, GT production system comprises of machining cells. 

2, Each component is classified and coded into family of 
components. 

3 * Each cell contains all machines and equipment needed 
to complete all the operations for the component families 
assigned to the cell. This is to reduce the inter cell 
movements of components. 

4, Each component has definite sequence of operations and 
the components move through the production cell with minimum 
back tracking. 

The design and implementation of a GT production system 
involves the consideration of the following problems. 

1. Component classification and coding schemes. 

2. Design of part families and machine cells. 
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5* Sequencing and scheduling of components belonging 
to part families assigned to the cell, 

4. Lot sizing of components. 

1.2.1 Component Classi fication an d Cod ings 

The objective of classification and coding schemes 
is to classify the conponents by their features so that 
components having similar code numbers possess similar fea- 
tijres* The three basic component features used for classi- 
fication are shape, function and manufacturing operations 
and tooling. Different classification schemes use different 
combinations of these features. Some of the important classi- 
fication and coding schemes are optiz system, VUOSCO system, 
PERA system and Brisch system (3). The problem with most of 
the schemes available in the literature is that they only 
consider the design aspects of the componont and do not account 
for certain important aspects like the machines used, 1die total 
requirement of component etc. Therefore, there is need to 
develop classification and coding schemes for GT production 
system based on features pertaining to production, design and 
resources used. Such classification and coding schemes would 
help in process planning, production control, data processing 
etc. of a GT production system. 

1.2.2 Design of Part Fa milies and Mac h ine, Cells? 

The design of groups and cells involves the identifi- 
cation of components to be made and machines to be insEfcalled 
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in each cell. Quantitative methods such as production flow 
analysis (PFA) or component flow analysis (CFA) can be used 
to establish the groups of components and machine cells. 
Malytical methods such as clustering technique, graph- 
theoretic approach can also be used to identify the compo- 
nent machine cells. These methods identify component machine 
groups by analysing the machine to machine routes followed by 
all Components. 

A layout for the machines in the cell has to establi- 
shed so as to madiine all the components in the part family 
with minimum back tracking. 

1.2,3 Sequ enc ing, and Scheduli ng of Part. Fam ili es; 

Having determined what is to be manufactured in each 
cell, the sequence of operations for the components must be 
decided. This may be arranged to achieve a particular aims 
maximum labour or machine utilisation, minimum throughput 
time, minimum setup time by sequencing components by simi- 
larity of tooling requirement, minimimi material consimiption, or 
an optimum combinations of these several objectives. 

Ideally, with GT the groups should receive a series of 
orders at regular period intervals. Under reasonable assump- 
tions such as all the components are processed in the cell, 
an extension of Johnson's algorithm can be used "^f or a two 
stage problem to minimise the make span time. Similarly for 
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a N-stage problem a branch and bound procedure can be used to 
minimise the make span time. The sequencing and scheduling 
decisions should incorporate the latest information in giving 
a solution. Also the supervisors should be given sufficient 
freedom to over-ride any solution as they can get the work 
done more efficiently. 

Loading on various machines and the cells should be 
uniform throughout the plant to avoid any misunderstanding 
among the management, workers and unions, 

1,2,4 Lot Sizi ng of Compo nents^ 

After determining the sequence of operations for each 
con^onent, then it is necessary to decide for every component 
how much to produce at each production facility. This quan- 
tity is known as lot size produced at the facility in a single 
s etup « 

1.3 MOTIVATION AND SCOPE FOR PRESENT STUDY s 

One of the important problems in the manufact\iring 
systems based on GT concepts is the determination of economic 
lot sizes for the various stages (machines) comprising the 
cell. The determination of lot sizes in GT manufaotiiring 
environment, is some what different from the traditional 
manufacturing system for the following reasons. In a GT 
production system, one has greater control over the movement 
of components compared to the traditional methods of 
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production (3). Fig. 1,2^ illustrates the differences in 
production throt:igh-put time for the traditional scheme of 
manufacturing and GT production system. The elements of 
throughput time in a functional production system are queueing 
time, machining time and transportation time. In GT produc- 
tion system the transportation time is completely eliminated 
because all machines needed for processing the entire part 
family are located within the cell. The queueing time is 
also minimum due to the greater autonomy given to the cell. 

So the throughput time for a given component of a part family 
mainly constitute the total time spent by the component in 
the GT cell. The total time spent by the component has the 
following two elements. 

1, Total setup time i.e. the setup time required on all 
the machines for the production of the given lot of components, 

2. Total machining time for processing the entire lot of " 
the component on the various machines in the cell. 

The total setx:^) time is constant for the component 
irrespective of the lot size. However, the total machining 
time of the lot would depend on the lot size and influence the 
throughput time for the lot. Pig. 1,3 shows the relationship 
between throughput time and lot size. The work-in-process (WIP) 
inventory increases with an increase in the throiighput time 
■vdiich in turn is a function of the lot size. If/IP inventory 




q ; queue I nq , ^ 

M : Total MQchini'rtq time 
T : Transportion time 



FIG. 1-2 DIFFERENGES IN PRODUCTION THROUGHPUT 
TIME OF FUNCTIONAE PRODUCTION AND GT. 







FIG. 1-3 RELATIONSHrP BETWEEN LOTSIZE AND 

THROUGHPUT TIME IN GT. PRODUCTION SYSTEM 
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cost is a major contributor in the total cost of production 
inventory system operating in a GT environment (4) , One of 
the ways to reduce the contribution of work-in-process inven- 
tory towards the total cost would be through the splitting of 
the lot to be manufactured at the machine/ stage into batches 
and transporting the batch to the next stage. With this 
particular idea in view several lot sizing models have been 
developed in the present thesis. 

The first model which is an extension of the basic 
Wilson's economic production quantity model has been developed to 
reflect the influence of characterised in cell type of 
manufacture. 

Since the GT production system permits the transporta- 
tion of the lot to the next stage in the cell hs no time and 
at lower costs, the lot can be split into batches while 
processing, and can be transported to the next stage, for 
further processing so as to reduce the work-in-process inventory. 
Under such an assumption, a constant lot size model with lot 
splitting is developed. A two stage algorithm, involving a ! 

heuristic procedure and an optimisation procedure, is presented. 

In a GT cell there can be madaines with high production 
rates involving high setup costs. To distribute the high 
costs associated with the seti;q), it may become necessary to [ 

' f 

produce lots of different sizes at various stages. For surh a 



13 


situation a variable lot sizing model assuming integet ratios 
of lot-sizes at successive stages is developed. The model 
is formulated using dynamic programming (DP) and the recur- 
rence relations are developed to solve the problem. 

Finally, the features of constant lot size model and 
variable lot size model are combined and another model called 
variable lot size model with lot splitting is developed. 

A he\jristic procedure is suggested to solve the problem. 

The solution methodologies of the various models have 
been coded in Fortran 10 for DEC-1090 system and were tested 
for randomly generated problems with input variations in a 
specified range. The computational experience with each model 
is reported. 

1,4 ORGMISATION OF THE THESIS" 

Chapter II deals with a brief literattire review on 
group technology with special reference to GT lot sizing. 

In Chapter III, four models viz., an extension of V7ilson 's 
EPQ model, a constant lot size model with lot splitting, a 
variable lot size model, and a variable lot size model with 
lot splitting, are presented. For every model, a brief 
statement of the problem is presented first and is followed 
by formulation and solution methodology. The solution 
methodology for constant lot size model with lot splitting, 
is explained using an illustrative example. 
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Further, computational experience based on solving a 
set of randomly generated problems of varied sizes is 
presented for the last 3 models. In Chapter IV, conclusions 
based on the present study along with suggestions for further 
work are presented. 



CHAPTER II 


LITERATURE SURVEY 


In this chapter a brief review of the literature on 
Group Technology with special reference to the lot sizing 
problem is presented, Waghodekar and Sahu (30) have given an 
excellent bibliography of more than 450 papers on the subject 
¥e present a review of the important literature on GT under 
the following categories. 

1. Design of cells and groups, 

2. Group scheduling and sequencing, 

3. Lot sizing problem, 

4. Performance evaluation of GT production systems. 

2.1 DESIGN OF CELLS AND GROUPS: 

The following approaches have been reported for the 
design of cells and groups. 

1, Rule of thiomb approach by Edwards (5) 

2, Composite component approach by Edv/ards (5) 

3, Classification and coding schemes by Burbidge (2), 

4, Flow analysis 

i) Production flow analysis by Burbidge (2), 

ii) Component flow analysis by El-Essaway (6), 
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5. Approaches using similarity co-efficients, 

i) Cluster Analysis by McAuley (7), 

ii) Graph theoretic approach by Ra;3agopalan and Batra (8), 

6. Cell formation using Monte Carlo simulation by Crookall 
and Baldwin (9), 

7. Mathematical classification by Purcheck (lO), 

8. Matrix clustering technique by Me Cormick (ll) . 

The first four approaches are some what qualitative 
in nature while the remaining are analytical approaches. 

A detailed discussion on these approaches is given by Wagho- 
dekar and Sahu (12). 

2.2 GROUP SCHEDULING AND SEQUENCING; 

Petrov (13) has developed four inter related scheduling 
models for different types of route sequences and component 
flows, Hitomi and Ham (l4) have suggested a technique for 
scheduling multi product multi-stage manufacturing systems 
using Ignall and Shrage branch and bound approach. Hitomi, Ham 
and Yoshida (15 ) considered group scheduling decisions under 
due date constraints. However all these models assume that 
the set-up time is included in the processing time. Kishore (l6) 
separated the setup time from processing time and developed an 
extension of Johnson’s algorithm for the two stage problem 
considering the criterion of minimising the make span time. 

For a N-stage problem, a branch and boimd procedure and a 
heuristic procedure has been developed. 
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In addition to group scheduling, machine loading and 
product mix decisions represent major problem areas for group 
production planning and scheduling, Hitomi and Ham (17) have 
considered problems from the view point of GT, for a single 
stage production, Agrawal (l8) has extended the work of 
Hitomi and Ham and suggested optimisation techniques for 
multi-stage problems. 

2.3 GT LOT SIZING', 

The lot sizing in GT is a special case of lot sizing 
in mxxLti-stage production inventory system. 

For a multi stage production system Crowston et al (l9) 
have suggested that in an optimal schedifLe the lot size at 
any given stage should be an integer miiltiple of the lot size 
at its immediate predecessors and suggested that the problem 
can be solved by examining all combinations of such integer 
values . 

Chakravarty (20) has considered the production planning 
and lot sizing problems , for mutually independaat machine 
component groups. Assuming the integrality theorem of 
Crowston (19) and no splitting of the lot for inter-stage 
shipment, production cycle time of every machine was found 
considering the set-i^) and inventory costs. Also a network 
based design approach to integrate the lot sizing and layout 
decisions has been presented. 
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Ignall and Veinott (2l) have suggested system myopic 
policies for multi-stage production system under continuous 
review with constant demand over infinite planning horizon. 
System myopic policies optimise a given obj active function 
^n.th respect to any two stages and ignore the multi-stage 
interaction effects. 

Wagner and Vttiitin (22) have developed a dynamic lot 
size model to solve the lot sizing problem of single product 
with known demand in discrete time periods. Zangwill (23) 
has suggested a network formulation for determining dynamic 
economic lot sizes with back logging. The network formulations 
facilitate the development of efficient dynamic programming 
algorithms for obtaining the optimal dynamic lot sizes. 

Goyal (24) developed a mathematical model for lot size 
scheduling on a single machine for stochastic demand, A 
method for computing the lot size in each time period is 
presented. Newson (25) developed a network based heuristics 
to solve the capacitated lot size problems with fixed resources 
and variable resources. 

2.4 performance MEASUREI'iENT BY SIMULATION; 

Gupta and Tompkins (26) have studied the performance 
of a GT production system with a simulation model written in 
Simscript, The performance characteristics include average 
stay time, intercell and intra-cell movements, number of 
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orders completed in time etc., Ang and V/illey (2?) have 
presented a simulation model which compares the Pure GT 
and hybrid GT» In hybrid GT, inter cell movements are 
permitted to certain extent while no inter-cell movements 
are permitted in pure GT production system. 



CHAPTER in 


LOT SIZING MODELS 


In this chapter the following lot sizing models have 
been developed, 

1. An extension of Wilson’s EPQ model, 

2. Constant lot size model vdth lot splitting, 

3. Variable lot size model, 

4. Variable lot size model with lot splitting. 

For every model, a brief statement of the problem is 
presented first and is followed by assvimptions, notations, 
formulation and solution methodology. The coaputational 
esqjerience for models 2, 3 and 4 based on solving a set of 
randomly generated problems is also presented. 

MODEL is 

3.1 an extension of WILSON'S EPQ MODELS 
3*1.1 Statement of the Problem s 

Consider a GT cell manufacturing a part family of 
components with known annual demand. The components are 
machined sequentially by the machines in the cell and the 
inter transfer time of components between machines is ne^igible. 
The problem is to determine the Economic Production Quantity 
(EPQ) for every component considering the costs due to work 
in process, seti^ and finished goods inventory. 
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3.1.2 i^sum pt ions,; 

1. All components are processed in the cell and inter 
cell movements are not permitted. 

2. Demand rate is constant over the time horizon. 

3. The inter transfer time of components between machines 
is negligible. 

3.1.3 NotatJ-PSi' 


R ; Cost of operating the cell/unit time. 

For a given component j of a part family, 

Q. ; Lot size, 

J 

D. ; Annual demand rate 

J 

M. : Unit material cost, 

J 

V. ; Value added in the cell, 

J 

W. ; Work in process inventorj^ value, 

J 

TS^ Total setup time for all machines used by the 
component, 

TM. : Total machining time on all machines used by 

O 

the component , 

h. ; Average inventory carrying cost. 

u 

3.1.4 Model_ Fc^ulatipn 


The throughput time per lot is given by.TGT = TS. + TM. .Q. . 

U u J 

As the component progresses through the cell, value will be 
added to it. The value added can be e3q)ressed as, 
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V = (Throughput time/unit) x Cost of operating the cell 

u 

TS . 

= (--qJ- -i- TM )R 

The work in process inventory value per cycle can be 
given as, 


VWIP/ cycle = (Unit Material Cost + value added) Lot Size 


The annual WIP inventory value can be written as, 


A. = ( VWIP/ cycle) x No, of cycles x Throughput time/lot 

= CM. (TSj Q^) 

Simplifying we get, 

A = D (M + i. V ) X (TS. TM. Q.) 

vJ vJO^vJ u vJw 

The finished goods inventory value per unit may be written as 

PG(Q.) = M. + V. 

0 o J 

The total cost is the sum of setup costs, finished goods 
inventory carrying costs and the V/IP inventory carrying cost. 

The total cost function can be written as, 

TC(Qj) = * hj FG(Qj) .. hjD.CM. + f ) 

(TS. + TM. Q.) 

0 0 <3 

Substituting for PG(Q.) and V we get, 

u J 
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TC(Q^) 


F.D. + h.D. TSt R/2 T 

Q.[ ^^h.(M. + TM. R) 
(at j 3 0 3 3 

RTM. h. TS. R 

+ h.D. TM. (M. + -'r . J--,.. J.,.. 

33 33 2^ 12 


+ h.D.M. TS. + h.D.R TS. TM. 
333 3 33 3 3 


(3.1) 


Since the function is a single variable convex function, 
differential calculus can be applied to solve for Q... Diffe- 

J 

rentiating with respect to and solving it by equating to 
zero, we get. 


=\ 


2F. D. + D. h. TS.^ R 

.... .A, ............. 

h (M. TM.R) + 2h.D.TM. (M. + (TM.R)/2) 
3 3333 3 


(3.2) 


Here Q.^' is the optimal lot size for the component 3 of the 
d 

part family. 


* 


Substituting <3. = Q. in Eq, (3.I) we get the optimal 
d d 

total CO st TC ( ot ) , 

d 

3.1,5 C ompar ison, with Basic Wilson's EPQ Mod^eli 


For the Wilson's EPQ model, the total cost expression 
without considering cost of the work in process inventory value 
can be written as 

F. D. Q. 

TCCQ.) = + FG (Q’J h. 

3 3 3 *^ 

3 

F.D. h. Q. TS. 

= [m 4- (— + TM.) R] 

Q. ^ a. ^ 

3 3 


(3.3) 



24 


The economic production quantity Q. can be obtained 

u 


similarly as, 


“a 


2F.. 

h. (M. + TM.R) 
0 1 0 


(3.4) 




Substituting in Eq. (3.3) the optimal cost TC(Q*.) can be 

u u 

obtained. 


The two models are compared for variety of problems 
with variations of input parameters within a given range. The 
ranges selected for the input parameters are given in Table 
3*1. The input data selected for a sample of 6 problems is 
given in Table 3.2. The economic production quantity and the 
total cost for the two models are presented in Table 3.3» Eor 
all the sample problems the consideration of WIP inventory 
results in smaller production lot size as well as total cost 
of the production inventory system, 

MODEL 2s 

3.2 CONST/iNT LOT SIZE T^TH LOT SPLITTING l 
3-2.1 Stat ement of ^ the_ Prob lem? 

Consider a GT cell comprising of N stages (each stage 
corresponds to a machine) manufacturing a part family of 
components with known annual demand. The components are 
produced sequentially on the various stages, in the cell. A 
constant lot is to be produced on all stages. The lot being 
produced at a particular stage can be split into batches which 
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Table 3.1s Input Ranges 


Variation in Variation Variation Variation in 

Demand in Material in Setup Machini:^ 

(pieces) Cost (Rs.) hrs. Time (minutes) 


1000-10000 2.0-10.0 1.0-12,0 8-35 


Table 3.2s Inpxjt Data 

R 75000 Rs./year, - 10 percent 

Production hours available for the cell/year = 2880 


Problem 

No, 

Demand 

Material 

Cost 

(Rs.) 

Total Setup 
Time TS . 

(Hrs.) 

Total Machini 
Time/Unit 

MS. (Minutes) 

cl 

1. . 

6000 

2.00 

9 

20 

2. 

4000 

4.00 

12 

25 

3. 

5000 

5.00 

8 

30 

4, ; 

LOOOO 

3.0 

11 

35 

5. 

9000 

2.0 

10 

23 

6* 

1000 

4.0 

6 

20 
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Total 3 m3’ Comparison of Two Models. 
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can be transported to the next stage for further processing, 
even -when processing of the remaining lot is still in progress 
at this particular stage. The inter transfer time between 
stages is negligible. The problem is to determine the cons- 
tant lot size for all the stages and the number of batches 
for the production of the lot at each stage such that the 
total cost of the system arising on account of setup, trans- 
portation and inventory is minimised. The lot size, the 
batch size at each stage and the number of- batches into which 
the lot is split must be integers. 

3.2.2 As_s imrptipns_? 

1, The demand rate is deterministic and constant over the 
planning horizon. 

2, For all the stages the setup costs are fixed, 

3, The inventory carrying costs are linear in nature. 

4, The transportation cost per batch at a stage is inde- 
pendent of the number of xmits transported. 

5, The inter transfer time of components between stages 
is negligible. 

6, The production rate at each stage is;: greater than the 
demand rate . 

3.2.3 Notat ion s 

For a component j of a part family, 

D ; Annual demand rate of the component (final product) 
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Q s Constant lot size, 

N ; Number of stages, 

For the component at any stage i, 
i Batch size, 

; Number of batches, 

; Setup cost per lot, 

h^ : Unit inventory holding cost per unit time, 

: Transportation cost per batch, 

: Machining time, 

: Elapsed time between stages i and i-i-1, 

R A real value R rounded to the nearest integer, 

R'- : A real value R roimded to the higher integer, 

R -j" 5 A real value R rounded to the lower integer, 

3.2,4 Mode l Fomulationj 

Fig. 3.1 represents the inventory building between 
stages i and i+l, for the case m^ > ^he first slanted 

line represents the emulative production at stage i. The 
comers of various triangles formed with this line indicate 
the availability of batches for transportation to the next 
stage i+l. The second slanted line represents the emulative 
production at stage i+l. The dotted lines crossing this line 
represent the depletion of stage i inventory. The trapezoid 
enclosed by solid lines represents the time weighted inven- 
tory at stage i* The inventory at stage i builds up over time 
period Qm^ during ^^hich y^ number of x^j^ sized batches are 
transported to stage i+1. 




FIG. 3-1 TIME-WEIGHTED INVENTORY AT STAGE i WHEN 
y =5; = 3 : AND m; >m|,-i • ...... 
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From the Fig. 3.1, we observe that for uninterrupted 
production, the elapsed time between stages i and i+1 is 
given by. 


The time-weighted inventory at stage i is given by the area 
of the trapezoid, 

A = I I- 

for 

Similarly when m^ ^iH-l» easy to verify thatj 

A = I [(x.m, + + x.m^ - Qm.))] 

The two expressions can be combined as, 

A = |- [2x^ min m^) + Ojm^ - “i+1^^ 

Substituting and simplifying, 


A = 


£ [ . 1 u n 

2 y,. 2 ‘“i ®i-!-l‘^ 


The average inventory cost can be calculated by multi- 
plying the time weighted inventory with inventory carrying 
cost .per unit time and the number of cycles. It is given 
by the following expression, 

N 

ifl ’’i 0 


Substituting for a, the expression for the average inventory 
cost is written as, 
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N min (m^ , m.) . 

D E h Q [ + i. im. 

i=l ^ ^ ^ 


m 


1+1 


13 


The total setup and transportation cost considering all the 
stages is given as. 


S(F^, T^ 


N F. 
D E ( 
i=l 


y.T. 

..1 + ll.A ) 
Q Q ^ 


The total cost function is obtained by summing up the setup 
costs, transportation cost and the average inventory costs. 
Thus the total cost of producing the component on the N~stages 
is given by, 


TC(Q,Y) 


K 

F. y.T, 

D E [1 

f 'll.. 1 

^ Q Q 

i=l 

+ I im^ 



+ h^Q 


, min 




where Y = {y^, y^} 

The problem of determining the optimum lot size can now be 
formulated as, 


Min TO (Q,y) = AQ + §. + I [( J-) y + -i- ] (3.5) 

sub; 3 ect to, 

y^ = Positive integer, 

Q = Positive integer 
Xf = ^/y± positive integer for all i 

A, B, a^, b^ > 0 


( 3 . 6 ) 
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where, 


A 

B 


N h. 

D S *5^- jm. - 
i=l ^ 


m 


i+1 


N 

D 2 F. 
i=l ^ 


and, 


b 


i 


*1 "lia (>»i. ^ i 

DT^, ¥ i 


5.2.5 Solut ipii jMethpdpl pgy ;. 

The optimisation of the objective ftinction given by (3.5) 
is carried out in two phases. In the first phase a heiiristic 
solution to the problem is found. In the heuristic procedure, 
initially the values of Q, and y^ are determined iteratively, 
maintaining integrality for each of them. In finding the 
optimal solution, the results of heuristic procedxn'e are lised 
as starting values for developing the upper and lower bounds 
on the lot size within which the optimal solution lies. Start- 
ing from the lower bound of the lot size the optimal solution 
is found by scanning several combinations of batch sizes and 
number of batches. An efficient scanning method is developed 
for this purpose. 

3.2.6 Heurijtic^ Prpcedi^^ 

The oboective function (3.5) can be rewritten as, 
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N N 

TC(Q,Y) = lA S (a./y.)]Q -J- [B 2 
i=l 1=1 

Substituting = (Vyj_ (3.5) we obtain, 


TC(Q, Y) 

= AQ + -A - + 2 a.x. -i- 

0 i=i ^ 1 

N 

i=l ^ ^ 


Let, 

ec.p = 

b. 

a.x^ + and 



0(7^) = 

(- 5 -’) Yi + (a^Q)/yj_ 



The total 

cost function given in 

( 3 . 8 ) can now 

be 

rewritten as, 




TC(Q, X) 

n N 

= AQ -f- + 2 0(xj 

i=l ^ 


(3.9) 

TC(Q, Y) 

= AQ -h 1 + ^2 0(y^) 


(3.10) 


i=l 


where, X = x^^} 

» 

Lemma iJ For P,q > 0 the positive integer K that minimises 
the function f(K) = PK + -^ is 

K = [ |- + 0.25]^/^ 2 

« 

Proof; The optimal K must satisfy, 

f(K - 1) f(K) (3.11) 


and, 


f(K) < f(K -r 1) 


( 3 . 12 ) 
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Using ( 3 . 11 )# we get, 

PK + < P(K - 1) + 

K “ K-1 

Rearranging, 

K(k-l) < q/P 

Adding 0«25 on both sides and simplifying, 

(k - 0,5)^^^ < f + 0.25 

(k - 0.5) < ( I + 0 . 25 )^/^ (3.13) 

Similarly using (3.12) we get, 

pk + - ^ < P(k+1) + 

K ^ K+1 

Rearranging, 

k (k + 1) > q/P 

Adding 0.25 on both sides and simplifying 
(k + 0.5)^ > (q/P + 0.25) 

(k + 0.5) > (q/P + 0.25)^/^ (3.14) 

Combining (3.13) and (3.14) we get, 

(k -- 0.5) < (q/P + 0.25)^^^^ < k + 0.5 

* 

This implies that K is rounded to nearest integer of 
(q/P -i- 0.25)^/^. 


This proves the lemma 
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Lemma 2; For P,q > 0 the minimum of f(K) = PK -!• q/K is 

K = (q/P)^/2 and 

f(K ) = 2(qP)^/^ 

Proof; Since f(K) is convex, for K > 0 solving *= 0, 

we get K and f(K ). 

For the time being relaxing the constraint on and 
using Lemma 1, we get from Eq. (3.9), 

Q = (B/A -1- 0.25)^/^ (3.15) 

Xf = (t>j_/a^ 0.25)^'^^ V i (3.16) 

Substituting Q = in (3.9) and again using Lemma 1 we get, 

(B/A) 0.25)^/^ ^ i (3.17) 

!)C • 

1 

After establishing X,Y vector from (3.7), a lot size \daich 
is nearest multiple of all s is obtained as, 

a = [(B -i- Z b y )/(A + E a./y )]-^/^ (3.18) 

i==l i=l 

Let L = Least common miiLtiple of Y-vector. 

q’ = (|;’*t) L 

With the updated value of Q the X and Y vectors are also 
updated as, 

x| = Q’/Yi, i 

and 

Yi = C + 0.25]^/^, V i (3.20) 

. X. ■■ . 
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With the available new Y-vector, the lot size Q ' and 
X-vector are modified using (3.18) and (3.19). With every 
iteration the total cost decreases and the heuristic proce*- 
dure stops when no further reduction in cost is possible or 
O’, X, Y vectors stablise. The various steps of the heuristic 
procedure are summarised below. 

Heuri stic Algorithm^ 

Step^ l; Calculate the co-efficients of A, B,aj^,b^. Set 1 = 1, 
TC = 38E + 11 (a high value). 

St^ep^,2,2 Calculate 

Q = (B/A 0.25)^/^ 'T- 

= (bj_/a^ + 0.25)^/^ V i 

y = [ .,1 (b/a) + 0.25]^/^ f V i 

^ X, 

1 

Step^J,: L = Least common multiple of Vector Y 

o' = [(B -I- E b.y.)/(A + E a./y. 

i=l ^ ^ i=l ^ ^ 

Q = ( L 

If (o' = 0) then Q’ = L 

^i = yI = [ -Tl- (B/A) + 0.25]^/^ t 

. N ^ , 

TC = AQ -i- B/Q' E 0(y. ) 

i=l ^ 
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Ste p 5» 


iS'S£. 


If q’ = Q and 

X;/ = Xf V i and 

Y;/ = y-f V i or 

TC' > TC GO TO STEP 6 
otherwise Step 5. 


Set I = I + 1 
Update Q =s Q , 


= y . V i 


X. = x' Vi 


TC = TC' 

GO TO STEP 3. 

Write Q, x^, y^^, TC as the heijristic solution, STOP, 


3.2,7 

The heuristic solution gives an upper bound on the cost 
of the optimal solution. Relaxing the integer constraint on 
in (3.9) s the minimum of 9(x^) can be given from Lemma 2 
as Then we have, from (3.9) 

TC = AQ + & + 2 E (a.b. )^/^ (3.21) 

^ i=l ^ ^ 

Substituting the heuristic solution cost as TC and 
sin 5 )lifj'-ing (3.20) we get. 


- 2aQ + p - 0 


( 3 . 22 ) 
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where , 

a = [ TC - 2 S 

i=l ^ ^ 

P = B/A 

By solving (3.22), an upper and a lower bound on the 
optimal lot size can be established as, 

Qy = a -!- (a^ - 

Qj = a - (a^ - 


Using as the starting point, the entire range of lot size 
is to be scanned for the optimal solution. 

For the development of an efficient scanning procedure, 
let us consider the effect of change in value of some over 
the given objective function. The part of the objective 
function influenced by y^ is given by, 

0 ( 7 ^) = (-^) Yi + (a.Q)/y^ 

Let us find the condition \jnder which a change in y^^ 
from y^ to y^ + 1 will result in 

0 (y^ + 1) <0 (yj_) 

Let us assume, 


0 (yj_ + 1) < 0 (yj_) 
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After rearranging we get. 


b, a.Q 

•S' i 


Thus, 


Q > [ ^- (Yi 




(3.23) 


From the above inequality, we infer that the starting 
points of ranges associated with changing to y^^ 1 for 

each stage separately are, 

qI = [ Yi^yi + (3.24) 

1 

In the scanning process, starting with Qj^ we first, 
establish the values for Y-vector using (3.16) and (3.17).. 

With the available Y-vector the lot size is obtained using 
(3.18). Then the new range of lot sizes for change in to 
y^^+l, for every i are established using (3.23). Tf the 
min (o'.) exceeds the upper bomd, we have reached an optimal 

jii 1 

solirtion, otherwise the y- corresponding to min (Q. ) can be 

1 ¥ i ^ 

changed to y^j^ + 1 as this decreases the value of 0(y^) , The 
entire scanning process is repeated again until all Q^' s fedl 
outside the upper bound The varioxis steps in the optimis- 

ing algorithm are given below. 
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St ep I t 




Set, 


a = 

P = 


TC = 

Q « 

Yi = 

[ TC - 

B/A 


^^^^heijristic 
^^^heuristic 
^heuristic 
**^1^ heuristic 

2 E (a b 

i=l ^ ^ 


¥ 

¥ 

2A 


i 

i 


= a - (a^ - 



%- 

a + (a^ - 


x’ = 

(b^/a^ + 0.25)^/^ t ¥ i 


y! = 

(*-f ) •:r V i 

X. 



1 


L = 

Least common multiple of Y— vector 


a' = 

N N -I/O 

[{(B + 2 b.y.)/(A+ Z a /b.)} + 0.25]-^/'^ 

i=l ^ ^ i=l ^ ^ 


Q‘ = 

i ^ i L 

Step 4i 

f 

Xi = 

a '/y'i w i 


TC ' = AQ' + b/q' -h Z c6(y'. ) 

i=l ' ^ 

If (TC*> TC) GO TO STEP 6. 
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Step. ^ : Set, TC = Tc’j Q = Q; 

^ % = x| V i 

Step_6: Oil = [ (y^ + V i 

M = min (Q"); j = K: Q" = M 

V i 1 ^ 

Step.2‘ TO STEP 8 

y ’ = y ' + 1; GO TO STEP 3. 

u J 

s tep 8 : Write TC, Q, y^, as the optimum results. STOP. 

3.2,8 Efficiency;, pi, the,. Scanni ng, P rocedt^^^^ 

The efficiency of the scanning procedure can be mea- 
sured in terms of the number of iterations involved as compared 
to the complete enumeration method. 

Let the feasible value of y . falls between y^^ and 

TT13X! 

y. . Let K. be the total number of feasible values of y-. 

.L aL* 

For a N-stage problem the number of iterations for con^lete 
enumeration would be 

le = 

In the scanning procedure used in the optimising algo- 
rithm only ‘one y^ is changed to y^+1 in every iteration and 
this can happen K.^~l times for each i. Thus for a N-stage 
problem, the number of iterations would be 

= i^j^-1) + (K 2 -I) + ... + (% - 1) -I- 1 

N 

= E K. - N + 1 
1=1 ^ 
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For N = 4; = 10 V i = 1,4 

I = 10^; = 37. 

Thus there is considerable reduction in number of iterations 
in the scanning process. 

3.2.9 Numeri cal E xamples 

A sample problem is solved numerically illustrating 
the various steps in the algorithm. The input data for the 
problem is given in Table 3.4. 


Table 3.4: 

Input data 

for the numerical example. 

D = 5000 

per year; 

Total 

hours available in a year •= 2880 

Stage 

m. 


h. T. 

1 

12.0 

18,0 

0.40 0.50 

2 

10.0 

15,0 

0.50 0.75 

3 

8.0 

16.0 

0.60 0.60 

Heuristic Solutioni 



Step 1: A = 

= 0.2566? 

B = 

245000; 

Pi = 

= {0.1157 

0.117 

0.1388}; b^ = {2500 3750 3000 } 


TC = 38E + 11 (a high value) 

S tep 2t a = 977; = {l47, 180, 147}? = {7, 5, 7} 

4 = 35? Q’ = 978? Q' = 980? x^ =*1140, 196, 140} 
TC' = 


618.1966 
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Step. 4i 618.1966 < 38E + 11; GO TO STEP 5. 

Step 3° 1=2; Q = 980; = [7 5 7}; = {l40 196 140} 

TC = 618.1966; GO TO STEP 3. 

Step^l; L = 35; Q' = 978; Q’ = 980; = {l40, 196, 140} 

Yi = {7, 5, 7}j TC’ = 618.1966 

Step 4; As Q' = Q; x'- = x. ; yl = y . and Tc'= TC the procedure 
stops. Write 


^ ^^heuristic 
(x ) 

^ heuristic 
Cv '■ 

'•^i^heuristic 

heuristic 


980 

{140, 196, 140} 
{7, 5, 7} 
618.1966 


Optimal Solution; 

S,tep^l; TC = 618.1966, Q = 980, y^ = {7, 5, 7}, 

x7_ = {140, 196, 140} 

Step^,2; a = 977.61218, p = 954793.45 
= 947.08, Qy = 1008.1429 

x^ = {147 180 147}, y^ = {6 5 6} 

Step^; L = 30, Q = 960 

Step,, 4; xj = {160 195 160} , TC' = 618.3785 
TC’ > fc, GO TO STEP 6 

Step_6i. = {952.638, 986.07, 952.7754}, M = 952.638, j = 1 
Step,_2- 952.638 < 1008,1429, y^ = 6 + 1 = 7, GO TO STEP 3 
Step.J.; Yj. = {7 5 6}, L = 210, Q' = 974, Q' = 1050 
Step 4; xj_ = {150 210 175}. TC' = 616.317, 

TC’ < TC , GO TO STEP 5. 
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TC' = 616.317, Q = 1050, = {150,210,175}, 

y\ = {7, 5, 6} 

Step. 6; = {llOO.Ol, 986.07, 952.7754}, a « 3, 

= 6 + 1 = 7, GO TO STEP 3. 

Step, 3.: L = 35, 0’ = 978, Q‘ = 980, TC’= 618,196, 
to' > TO, GO TO STEP 6. 

Step,^; = { 1100. 01, 986.07, 1100.17}, 3=2, 

= 5+1 = 6, GO TO STEP 3. 

S tep 3 -. O' = 990, Q' = 1008. 

Step,4j 1^^}» TC'= 618.269, 

TC’ > TC, GO TO STEP 6. 

St.ep.,.,6i Q^' = {1100.01, 1166.7587, 1100.17}, M = 1100.01, 

M > Q^, GO TO STEP 8. 

Stjp_ 8 ; Writ e , 

(«opt =1°50, 

= *^50, 210, 175}, 

(yi^opt = t7, 5, 6} 

3»2,10 Co mputa tio nal Exp erience: 

The algorithm has been codad in Fortran~10 and imple- 
mented on DEC-1090 system. Number of problems of varied sizes 
(Number of stages) were tested. The input parameters viz., 
demand, setup costs, inventory holding costs and transporta- 
tion cost, were selected randomly. The ranges selected for the 
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various input parameters .are given in Table 3.5. For each 
Table 3.5! Ranges of input data. 


No, of 
stages 


m. 



T 


i 


1-10 1000-15000 5-40 10-40 0.10-0.80 0,25-3.0 


problem size ten problems were solved. It was observed that 
in most of the cases the heuristic solution was obtained in 
less than three iterations and in no case it exceeded five 
iterations. In all the cases except two cases, the heuristic 
solution was fo-und to be optimal. One of the cases in which 
the heuristic solution was not found to be optimal is given 
as an illustrative example in the previous section. 

Though the heuristic procedure gave the optimal solu- 
tion in most of the cases, the optimality could not be 
guaranteed. Further, it was found that of the total CPU 
time for solving a given size problem, the heuristic procedure 
consumed less than 50 percent of the time. It optimality 
is not be guaranteed, considerable saving in the computa- 
tional effort may result in by simply using the heuristic 
procedure to obtain the solution of the problem. 

The effect of number of stages on the compiitational 
time was investigated and is presented in Table 3*6. For 
single and two stage problems the con5)utational time was found 
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Table 3 •St Con^jutatlonal performance. 


Ko. of"* Averai^^CPtF time'*" 

s tages i n milli sec . 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


92.5 

75.0 
34.2 

36.5 

40.1 
53.25 
74.8 

154.4 

226.5 
645.25 


to be more than that of a ^-stage problem. This can be 
ejqplained because in the single and two stage problems, the 
number of combinations of X and Y vectors to be evaluated 
is greater than those for the 3-stage problem, fbwever, 
for problem of size greater than 3-stages, the computational 
time is found to increase with number of stages. This is 
due to higher number of enimierations to be carried out to 
encompass the total n\jmber of stages in the problem. This is 
shown graphically in Fig. 3.2. . 

MODEL 5 

3.3 VARIABLE LOT SIZE MODEL: 

3.3.1 State ment o f the Problem : 

Consider a GT production cell comprising of multi- 
stages where in the lot size at each stage is an integer 



CPU Time in milli seconds 



FIG. 3-2 COMPUTATIONAL PERFORMANCE OF CONSTANT 
LOTSIZE MODEL WITH LOT. SPLITTING 
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multiple of the lot size at its succeeding stage. The 
demand rate is constant over the planning horizon. The 
problem is to determine the optimal lot size for each stage, 
so as to minimise the total cost of the production/inventory 
system. 

3 *3 *2 ^ 

1, The production rate at each stage is greater than the 
demand rate , 

2. The demand rate is constant over planning horizon. 

3* The setup costs are fixed at each stage. 

4, The lot is transported to the next stage onlj'’ when 
the entire lot is processed, 

5. No shortages are allowed. 

3.3*3 Not ati on; 

D ; Demand rate of the component of the part family 
N ; Number of stages 
At any stage i 
Q. ! Lot size 

I Fixed seti^ costs 
hj^ ; Inventory holding cost per xmiit 
m . ! Machining time 

wim 

i Transportation cost per lot 
: Positive integer >, 1 
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• Cost function with lot size = %^i^ 
tj^(Ki,Qjj) ? Transfer function of and 

3* 3*4 Models j Forp iulatipn: 

For a multi-stage production system the inventory at 
a stage i is defined as the number of units which have passed 
through the stage i but not left the system. For such a 
system Crowston et al (19) have shown that the lot size at each 
stage should be an integer multiple of the lot size at its 
succeeding stage. This is known as integrality theorem. The 
theorem is stated below, without proof. 


Integrality Iheorem^ If i denotes any stage, a(i) denotes its 

successor stage and N denotes the final stage then there exists 

a set of optimal lot sizes Q 2 , ..., such that for all 

i < N the ratios, 

Q. 

K. = « are positive integers. 

^ \(i) 

The proof is given in Crowston et al (19). 


The cost f motion fCCj.?, K^) at each stage consists of 
setup costs, inventory holding costs and transportation cost. 
For each cycle, the setup and transportation cost for stage i 
will be (Fj^ 4- The total seti^ and transportation cost 

of the stage for fulfilling the demand D in lots of will be 
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Fig. 3‘.3 represents the inventory huildij^) at stage i. 
The shaded area represents the time weighted inventory at 
stage i. The shaded area is given by 

A *= I 0^ ( 

The average inventory holding cost is obtained by multi- 
plying the time weighted inventory vdLth number of cycles and 
the unit inventory holding cost. Average inventory carrying 
cost is given by the expression. 


' I h. Q, (l-Dmp 

The total cost function for stage i is given by, 

D(F. -i- T.) 

f(%, K^) = — + I h. (l-Qi. ) 

D(Fi + T.) 1 

= (l-Qnp 

^i % 

Ihe objective is to minimize the total cost for all 
stages. This can be written as, 


TC = Min {f(%, %) + f(Qij, + ...+ f(% %)} 
s/t = K^Qjj for i = 1,2,.,,, M-1 (3.25) 







The problem given by (3.25) can be solved using Dynamic 
Programming (DP), D 3 mamic programming can be used only if 
the cost function is decomposable, mathematically, the 
total cost function is decomposable if it satisfies the follow- 
ing theorem. 


Decomposition Theorem: If a real-valued return function 
0j^( f ^ 2 * * * * • satisfies 

a) separability condition, i.e,, 

^2* **** ~ ^N-1 

where is a real valued function. 

b) 01 ^ is monotonic cro'iSfiiiTig function of 0 j^^ 2 . 

every fj^, then 0^^ is said to be decomposable. 

Since the total cost function given in (3.25) is separa- 
ble and monotonically non-decreasing, it is decomposable. 


The d^mamic programming model formulated can be re- 
presented diagramatically as shown in Fig, 3.4, The recursive 
relations can be written as. 


Ki-i 


= (K^, 


(3.26) 


TC^ = f’in ^ (T^i) + f(Q„. i%) (3.27) 

all prede*- 
cessors 

^^N-l = ^%-l^ (3.28) 

all pre- 
decessors 

I t 

r" -5 
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From these recurrence relations, the total cost fxmc- 
tion can be optimised, using forward recursion. It involves 
optimising the final stage given the initial condition that 
Kjj = 1, Then the proceeding stage is optimised for 
Qj^. For the available K-vector an optimal lot size Q.-^ is 
fomd. This continues till the final stage is reached, 

3'* 3.6 Computation al Exper ie nce ; 

The dynamic programming algorithm has been coded in 
Fortran -10 and implemented on DEC-10 system. Number of 
problems of varied sizes (number of stages) were tested. The 
input parameters viz., demand, setup costs, inventory holding 
costs and transportation cost were selected randomly. The 
ranges selected for the various input parameters are given 
in Table 3.7. 

Table 3.7^ Ranges of input data. 


No. of 
stages 


D 









1 - 10 1000-20000 5-40 10-40 O.IO-O.SO 0.25-3.0 

The average computational time required for solving 
problems of varied sizes (in terms of nimiber of stages) was 
investigated. For each problem size, five problems were solved. 
Table 3.8 gives the average computational time requirements 
which are presented graphically in Fig. 5.5. It is observed 
that the computational time requirements vary ejispoiicnttally 
with the number of stages. 



CPU Time in milli seconds 



FIG. 3-5 COMPUTATIONAL PERFORMANCe OF OP 
ALGORITHM 
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Table 3.8s, Average computational time vs. number of 
stages. 

No, of stages Average computational 

time in milli sec. 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


18.9 

29.1 

34.2 

45.0 

54.1 

66.4 

80.3 

96.5 
119.1 
144.6 


MODEL 4s 

3.4 VARIABLE LOT SIZE MODEL WITH LOT SPLITTINGS 
3.4.1 St at ement_ of, the. Problems 

Consider a GT production cell where in the lot size 
at each stage should be an integer mifLtiple of the lot size 
at its succeeding stage. In addition, the production lot 
at each stage can be split into batches. A batch can be 
transported to the next stage even when the lot to which the 
batch belongs is still being processed at the current stage. 
The lot size, batch size and number of batches at each stage 
should be integers. The problem is to find the variable lot 
sizes, the number of batches and the batch size at each stage 
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such that the total cost of the pro duct ion-inventor^'' system 
for the cell is minimised. 

3.4.2 As s mipt i pns » 

1. Demand rate is constant over time horizon. 

2. Inventory carr^ring costs are linear in nature. 

3. Inter-transfer time is negligible. 

4. The lot size at stage i is restricted by the capacity 
of the stage. 

5. The maximum batch size depends upon the load carrying 
capacity of transport equipment at that stage. 

3.4.3 Not atipn^c 

For a given component of a part family, 

D s Demand rate, 

N s Number of stages. 

At any stage i, 

• Lot size produced, 
m^^ s Machining time, 

h^ i Inventory carrying cost per unit, 

t Batch size, 

y. s Number of batches, 

•'1 

gj^ ; Maximum load carrying capacit 3 r, 

: Maximum lot size permi'fcted, 
s Integer ratio, = ^i-l/%^ 

Ei_^i: Earliest time at which production can be 
started at stage i-fl. 
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: Real value R rounded to higher integer, 

Rsl/i Real value R rounded to lov/er integer, 

R$: Real value R rounded to nearest integer. 

3'. 4 . 4 del. For mul at i on: 

Fig, 3.8 shov7s the inventory b\aildup at stages i and 
i-!-l, when the machining time at stage i is greater than at 
stage i+l. The iJipper slanted line represents the uninter- 
rupted production at stage i with a slope of l/m^. From this 
stage number of equal sized batches 
transported to the next stage i+1 as and when they are 
completed at stage i. 

Since m^^ > production at stage i-:-! cannot be 

started as soon as the first batch arrives at stage i-i-l. There 
should be some elapsed time after which only production at 
stage i+1 can be started. The second step function represents 
production at stage i+l, Hov/ever, the second slanted line 
V7hich appears below the first one, should satisfy the conti- 
nuous cumulative demand. The minimum earliest start time for 
stage i-Hl woiiLd depend on \i7here the tv7o step functions meet. 
This is necessary to satisfy the condition that there can not 
be production of the component at stage i-!-l without its being 
processed at the earlier stage i. Using this condition, the 
total elapsed time for both the stages can be obtained. Let 
j, 0 = {l,2, .. ., y^. } represent the batch at stage i v/hich is 



L0T5IZE 



FIG. 3-6 TIME-WEIGHTED INVENTORY BETWEEN STAGES 
i AND i + 1 WHEN 
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currently being processed at stage i+l. The total elapsed 
time at stage i is given by, 


ET. = (i-i-l) 


CiBii 


(3.29) 


The total elapsed time at stage i+1 can be e^qjressed as, 

1+1 1+1 


ix.m. 


(3.30) 


Equating (3.29) and (3.30) and simplifying, the earliest start 
time at sta.ge i+1 is given by. 


®i+l 


0 ^ ^ 3. 

(o+l) x.m. - 

The condition that there can be production at stage i+1 only 
v;hen stage i supplies it implies that the 

elapsed time at i+1 >_ elapsed time at i 

j !/C « 

=i? * 3-1 arj i- Qi.,1 (J- - (j-a) Xj^m^ 

1 “r_L 

for 0 < j yj_~l 

Rearranging and simplifying we get. 


E. = x.m. 
1+1 1 1 


ma3c 


o<d.<y-i-i 


L0(d)J 


(3.31) 


where , 


0(d) = dx^ ^ ^i+1 "“i+l 




For the case v;hen m^^ production at stage i+1 can 

be started as soon as first batch comes out from stage i. 
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This is because machining time at stage i is less compared 
to that at stage i-i-1 and so there will be enough units at stage 
i+1 for continuous prodioction. 

Hence the maximum earliest start time would be, 

\-i-l " ^i®i 

Since (mj_ ~ is alv/ays non positive 0(j) equals zero 

hence the expression (3.3l) can be used to determine 

To find the average inventory holding cost, let us 
first determine the time-weighted inventory. This is given by 
the shaded area in the Fig, 3,7. The area of the shaded 
portion is found by subtracting the areas of triangles formed 
from the trapezoid. 


Area of the trapezoid is given by, 
A = 


^i 




Area of the triangles can be v/ritten as, 

^ ^i-:-l ^i-i-1 ^ fe “i^ Q 

1t± 

Therefore, the shaded area can be expressed by, 

I % m^j - 


M 


MiJltiplying the time-weighted inventory by the inventory 
holding costs h^ and number of cycles, we get the average 
inventory holding cost expression as, 

I % ^^^i-i-l J-- mj.) - iii 
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By adding the total setup costs and total transportation 
costs, the total cost e^q^ression can be v/ritten as, 

1 


K F. b.T. 

TC = D 2 [( .1 ) -, 1 h. [2E._,, 

i=l i 


^^i ^ D- “i^ - ^i-M ^ IS’ - ^i+l^'>^^ 


(3.32) 


¥e have the follov/ing constraints to be satisfied at 
stage i, 

a) The batch size at stage i cannot exceed the maxiiiiuia 
load carr 3 ?-ing capacity at stage i, 

for all i = 1, , . . , N 


X,* g • 


b) 


The lot sizes follov/ the Crov/ston integral! 'by theorem (l9), 


Q 


i-1 

A' 


for all i - 2, , , . , N 


c) The lot size at stage i cannot exceed the maximum 
lot size permitted at stage i, 

Qi <_ for all i = 1,..., N 

d) The lot size is to be divided into equal sized batches, 

^i “ *^i/^^i* integer for all i = 1,..., N 

Using h^ = 0, the complete optimisation problem can be written 


in terms of as given belov/: 


N P 

liinimiSG TC = D 2 i. 'S' — m^^y (hj^’“h^^_2^j/2 

i=l i 


T. 

h. E. - + ■ A 
1 i-:-! x^ 


(3.33) 
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Hence for the above value of n, E. - would be minimum. 

This gives the lov^-er bound on E. Substituting 

Q. 


®i+l 


t - 1 in (3.31) v/e get; 

Q. 

x.m. 0 (.-y*A-|x_ 1) 


Lemma 3 5 For a real value R, we have, 

f ) 1 - 0 

I R4. ) 4 ^ - 0 

Proof: Always would be less than unity. Rounding to 

the lower integer always results in zero. Using the 
above lemma, and simplifj^ing the expression for 


we get; 


1) (ffli - ai+i)] 


(3.38) 


The two lov^er bounds for both cases can be combined 


into one b3'’ introducing a variable as follows-: 




(3.39) 


where , 


i 1 if m^ < 


I " kV" ®i - ®i+l 


3, 4, 5. 2 Bounds on Total Cost: To find a lower bound on total 
cost, the integrality restrictions on S^, y^ are relaxed 
and lower bounds on the earliest start time are substituted 


for E, . 
i*i’i 
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Relaxing integer restriction on the expression 
for is given by, 

'i> -A -'•5‘i •" '“iti 

(3.40) 


v/here, 


i 0 ®1 “i+l 


i " ■ 1 if 

Substituting (3.40) in (3.31) and writing all the terms in 
terms of we get lower bound on total cost TC 

N F, . 

TC' = ° t qY + Qj, {(^- - m^:,(h,-h^)/2 

*'i-l ^“i-1 - ■" "=i^i 

T. 

^ ? 4 m . T } + J 

JL .1. jK. * 

1 

This can be wit tan as, 

N A. G. 

TC = D f ' aj. * B.Q^ * f ] 

where 


^i 

= F. 

1 


B. 

1 

- L(fe 

- ”i)^ Wi 

TT 

= h.L( 


^i 

= T. 

1 



W‘'“i'' 


Hence the relaxed version of the original problem can be 
Witten as. 
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I'linimise fc 

(3.41) 


lower bound 
feasible 
by the heuris- 


in the previous 
s. With this 

the best integer ratios s are established, V/ith the given 
values the new lot size at all stages are modified. The 
number of batches at each sta.ge, are established using Q^, 
values, by an efficient search method. 

With a complete solution available i.e. S, Q, Y vectors, 
the lot sizes are modified using the eq. (3.32). The entire 
procedure is repeated until no further reduction in cost is 
possible. The algorithm converge as the total cost function 
is a non^decrsasing function and there are finite combinations 
of S and Y vectors. 



^i 

V 

i 

— !#•••> 

N 

Xi < 

Si 

¥ 

i 


N 

a, < 

% 

¥ 

i 

If ♦ • • f 

N 

^i < 

^i 

¥ 

i 


K 

Qf, x^ 

> 0 






The above problem can be solved to find a 
on the cost of the original problem. Using this 
solution close to the lower bound is established 
tic procedure described in the next section. 

3.4,6 Heuristic Procedure; 

3.4, 6,1 Outline of the procedinre; 

The relaxed version of the problem solved 
section gives a lower bound on cost and a set of 
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3.4«6,2 Development of the Procedure* 

Initially, vvith the available solution for the relaxed 
problem, the values can be established by rounding off 
procedure, 


~ L 


‘^N^'\+l \-:-2 %I+1-'' 


(3.42) 


for K = N, N-1, 3,2 


given 


= 1 


The maximm allowable Qi^ due to the lot size restriction 


at each stage is given by, 


So the permissible lot size at Final stage ^ 


(3.43) 


Qjj = min (Q^jt, Qj^; 

Given and S-vector the nev/ range of Qj|_' s are given by, 

'^’^K ^K-J'l ’** ^N-i-l^ for K s= 2,3,«»«, J'<“1,N 

(3.44) 


To find the Y-vector consider the following cases* 
a) Ylhen = ^i-i-1 

In the total cost eq. given b^?- ( 3 . 32 ), convert each 
X. into Q./y. and by v/riting the total cost in terms of 
y.. , we get, 

T. Q. 


(3.45) 



"When m. < 

m. T 

C. - 0 

X ' ' 

l+l 



T. 

h.Q, 

e{y,_) = 

(a- 

y . + • — 

yi 

s • is • 

vi ; 



m. 


Since Q (y^ ) is convex function using Lemma 1 in section 3.2, 

Q 2 

n = E H. 0.25]''/2 t 


(3,46) 


To satisfy the constraint, we have, 





. *^i ■ 


^i 


max vy^ , — T 

Case 

(ii) 

I'ftien m^. > m_. v' _. 



X ±->±j ^ X 


Using the Lemma 1 in s; 

lines 

that 




^i 

rr 

L * ■ ’ hj^ ^i-’i **' 

and 



0.- 


^i 

c=: 

max [y., 

^ ®i 

h) 

If 


'i “ ^i-'-l 


(3.47) 


1^/2 ± (3.48) 


(3.49) 


Othervrise the O(yj^) would be modified from (3.32) as given below 


T. 

1' 


©(y,* ) = (q 0 y.- -5- K- 

^ Oi 1 1 




gx. 

a/r4' (1/0 - 

1.-M 


(3.50) 


Let y^ be the value of y^ when 
already noted that 


we have 
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Therefore, from Eq. (3.43) wo have* 


T. 


h . Q . m . 


> (q^-) y-i ^ ^ 


‘ Q ^' ■'1 




(3.51) 


By solving inequality given by (3.5l)j, two roots of y^. can 
be obtained. Searching between the two roots for minimum 
©(Yj^) would give the optimum value of y^. 

Nov/ one complete set of solution viz., 

Cl--Vector, S~vcctor and Y-vector is now available. This 
can be substituted in Eq. (3.52) to get a better value of Qj^. 

Simplifying (3.32), we get^ 


TC = 


V/ 


a. 


ZQ, 


where, 


N 


N 


N 


(3.52) 


« = D (Fi ■> yiq)/(s„,s,,,^^ 

Z = D ,5; lSjj ... (bj_ - 

1=1 

(l“i- ) 

■'■^i-l ^i-^l ^“i~l “ ^ 

oi^\*i-l 'T 


By solving the e^q^ression given in (3.45), we get. 


Q 


(| V 0.25)^^:^ 


(3.53) 


The iterative process continues till the lot size obtained 
in (3.53) stabilises. 
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The steps in the algorithm can be smmarised as belov/i 

Solve the relaxed constraint problem given in (3«4l). 
Establish the values from the eq,(3.42) . Modify the 
lot sizes according to equations given by (3.43) and 
(3.44), I = 0. 

,Step.,3/" I » I + 1 , if m.. < m. ,T, calculate y. -value 

according to equations (3.46) and (5.47), otherwise, 
calculate - values according to equations (3.48) 
and (3.49). 

St ep^ 4 1 If = ^i+1’ Step 3» otherwise calculate the 

roots of y^ according to inequality (3. 51) and search 
for the optimum y^ between the two roots, GO TO STEP 3* 
Step^ J_; Calculate Q.-^ using eq. (3.52). 

Step_ 6.; If Qj^ = Q, Go to Step 7» otherv/ise = Qj.ji ^ 

GO TO STEP 2. 

Step„7^: Write the heuristic results, stop. 


3.4.7 Cpmput atipnal, E 2 q>e ri eiic e> 


The variable lot size model with lot splitting has been 
coded in Fori:ran-10 for DEC 1090 system. Number of problems of 
varied sizes (Number of stages) were tested. Kie input para- 
meters viz,, demand, setup costs, inventoiy holding costs and 
transportation costs were selected randomly. The raises selected 
for various input parameters are given in Table 3.9. For each 
problem size five problems were solved. It was observed that 
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■the heiiristic 301111:1031 was obtained in less than five itera- 
tions and in no case it exceeded 10 iterations. 

The effect of number of stages on the computational time 
was investigated. Table 3.10 gives the average computational 
time requirements which are presented graphically in Fig. 3.9. 


Table 3. 9* Ranges of input data 



1 - 10 1000-60000 2-40 10.0-40.0 0.1-0. 8 O.25*'3.0 


Table 3.10: 

Computational performance 

No. of stages 

Average CPU time 
milli seconds 

1 

20.1 

2 

31.5 

3 

39.5 

4 

56.0 

5 

70.0 . 

6 

85,1 

7 

107.3 

8 

129.8 

9 

161.7 

10 

185.3 


3.4.8 Comparison, of Jfodel^ 

The performances of Ilodel 4 and Model 2 have been com- 
pared for problems of varied sizes. For each problem the 




RG.3-9 COMPUTATtONAL' PERFORMANCE OF VARIABLE 
LOT5IZE MODEL' WITH LOT SPLITTING 


73 


randomly generated iE 5 )ut data was kept the same. The compa- 
rison was based on the total cost of operating the production/ 
inventory system and the computational time. From Table 3.11 
we observe that the total cost of the inventory production 
system is lower for the variable lot size model with lot 
splitting (Model 4) as compared to the constant lot-size model 
with lot splitting. However, the results v/ere found to be 
otherwise for single stage problem. This might have occured 
due to the heuristic solution procedvire followed for model 4, 

On an average, based on the problems considered, a reduction in 
total cost of about 10 percent was observed for the Model 4, 

Since the amount of reduction would be input data depen- 
dent, the only conclusion one can drav/ is that model 4 shoifLd 
be preferred over model 2, Further, model 4 requires lesser 
computational time for the same size problem as compared to 
model 2 as is evident from Tables (3.6) and (3.10). 
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Table 3,ii 


Comparison between Models II and IV 


Problem 
Size 
(No. of 
Stages) 



Percentage P^eduction 
in cost for Model IV 


705.1596 700.029 

1311.1646 1409.087 

1190.2946 1339.683 

1596.0257 1786.867 

1397.6391 1724.614 

2086.1632 2184.217 

2078.8431 2420.355 

2365.1158 2772.907 

3426.9020 3781.078 

3433.4674 3901. 995 


“ 0.7342 

6.949 

11.151 
10.6802 
18. 96 
4.489 
14.11 
14.70 
9.367 

12.007 



CHAPTER IV 


CONCLUSIONS AlO SCOPE FOR FURTFIER STUDY 

4.1 CONCLUSIONS? 

In this thesis, we have developed, mathematioal models 
and solution methodologies for lot sizing in GT production 
system. The following conclusions can be made on the models 
developed and the solution methodologies. 

1, The consideration of WEP inventory in GT production 
system results in lower lot size and lower total cost for the 
production of a given component of apart family in the cell, 

2, The splitting of the lot into batciies reduces the t/IP 
inventory significantly. 

3. The performance of the heuristic procedure in model 2 
is quite encouraging in terms of compirfcational time and its 
abi].ity to generate optimal solutions. 

4. The variable lot size model with lot splitting v/ould be 
preferable to the constant lot size modal with lot splitting. 

4.2 SCOPE FOR FURTHER STUDY. 

The models presented in this v/ork 4ave been developed 
based on several assumptions. These assuugotions can be relaxed 
to make the models more realistic. Specific models need to be 
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developed for the following situations. 

1. Demand instead of being constant varies with time, 

2. Demand is stochastic. 

3. Stages comprise of more than one machine and the 
capacity of the machine at each stage is limited. 

U Setup costs are sequence dependent. 
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